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Abstract 

The proper-time 4d path integral is used as a starting point to 
derive the new explicit parametric form of the quark-antiquark Green's 
function in gluonic and QED fields, entering as a common Wilson loop. 
The subsequent vacuum averaging of the latter allows to derive the 
instantaneous Hamiltonian. The explicit form and solutions are given 
in the case of the qq mesons in magnetic field. 



The path-integral formalism in quantum mechanics, created by Feynman [TJ 
[2] is an important benchmark in the development and in our understanding 
of quantum theory. Many varieties of this formalism and new methods to 
solve the problems, which seemed before unsolvable, have been suggested 
since then, see the books and review papers [31 HI El El [TJ El [9] , summarizing 
the modern achievements in this field. 

The extension of the path-integral formalism to the quantum field theory 
was done in several directions. One of the most known line of development 
was started already in [2], where field variables, i.e. electromagnetic poten- 
tials Aj(x, t), (p(jc, t) play the role of quantum spacial coordinates q(t), and, 
the resulting path integral is becoming the functional integral. This line is 
now a part of standard lore, present in many textbooks, see e.g. El H El E] . 

It is however important, that in all works of this direction the path inte- 
gration concerns spacial coordinates and/or field variables, but not the time 
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coordinate, and in this way one can say, that this development is similar to 
the path integrals in quantum mechanics, where time plays the ordering role 
and stays outside of the realm of fluctuating variables. 

Another and more general approach, unifying space and time coordinates 
in the path integral, is based on the proper time coordinate. The latter 
was introduced by V.Fock [TO] and J.Schwinger [TTJ, who used proper time 
formalism for the field theory in external electromagnetic fields, however did 
not exploit path integrals. 

In QED path integrals, based on the proper time, were suggested in [12] 
and developed in [T3] . 

Path integrals for QCD both in time and space variables, using proper 
time as an ordering variable were suggested in [HI [T5J [TH] . The first use of the 
QCD path integral for quarks and gluons was done in [17] and exploited to 
demonstrate the confinement due to field correlators (stochastic confinement) 
for a review see [18] . 

The full form of the path integral in QCD for quarks and gluons, based 
on the proper time ordering, was given in [T7J for T = and in [T^] for 
T > 0, and different approximations were reviewed in [20] for some relativis- 
ts models and in [21] for QCD. It was called the Fock-Feynman-Schwinger 
representation (FFSR) and we retain this name in what follows. 

Based on FFSR a new relativistic Hamiltonian was derived in [221 123] 
for quarks and in [25] for gluons, where a new important variable was 
introduced - u, playing the role of the einbein variable [2E]. Its average value 
ujq is the average quark (or gluon) energy and explains the appearance of the 
notion of constituent mass in earlier models. The relativistic Hamiltonian 
with einbeins Ui allows to calculate all low-lying states in QCD: mesons, 
glueballs, baryons and hybrids from the first-principle input: current quark 
masses, a s and string tension a, see [23 ESI [23 for reviews. 

However, the introduction of Ui as einbein variables, being successful, 
is an approximate procedure, and its limitations and corrections were not 
enough clarified in the literature. An attempt in this direction was done in 
[3D] , where the fluctuation of the time coordinate in the path integration was 
substituted by the fluctuations integration in Aw,. The resulting expressions 
for quark decay constants of mesons in [30J are quite successful in comparison 
with experiment, however the exact scheme of approximations was not clearly 
stated. 

An additional impulse for a development in this area was given recently 
by the inclusion of high magnetic field B in the dynamics of QCD and QED, 
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see [3TJ E2J GSJ [Ml ES] for a recent papers. In this case U{ depend on B and 
might vanish or grow fast (depending on quark spin projection), which calls 
for a careful analysis of all corrections. 

In the present paper we undertake such an analysis and rederive different 
forms of path integrals and relativistic Hamiltonians for the QCD and QED 
systems, typically for quark-antiquark or atoms, taking into account both 
QCD and QED dynamics in the first case. 

The first thing we meet confronting Ad path integrals, is the problem of the 
time-coordinate fluctuation, which necessarily requires distinguishing average 
(ordering) time and fluctuating time, similar to the old notion of the Zitter- 
bewegung. We analyze this phenomenon, comparing the Bethe-Salpeter and 
path-integral formalisms and show how the latter can be developed using the 
fact, that all dynamics is contained in the Wilson loop formalism augmented 
by spin insertions. 



2 Path integral: treating time fluctuations 

We start with the simplest example of a scalar particle in external field, this 
problem was considered for QED by Feynman in [12]. 

The scalar one particle Green's function is (in Euclidean space-time) 

y) = ( m^-D^ I = r ds ( D * z) *y e M~K)Hx, y) (1) 



A 1 / xy 



where = - ieA^ 



$(x,y) = expie / A^dz^ (3) 

Jv 



and 



(D A z) - lim TT l 3 -^- [^A^n=An)-t*-v)) Ne = s U) 
1 )xy ~ nZz^J (ATiefJ (2vr) 4 ' 1 ' 

At this point it is important to stress the difference between the nonrel- 
ativistic quantum-mechanical and relativistic path integration: in the first 
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case one has (D 3 z) = (D 3 z(t)) in (pQ) and the time variable t has the ordering 
character: the consecutive pieces of trajectory z(t) are ordered by time. In 
the relativistic path integral this role is given to the proper time r, s while 
the "time" z 4 (r) is fluctuating together with spacial coordinates z(r). In 
terms of any local field theory and Bethe-Salpeter type of equation this is 
allowable and necessary, since any moment of time z 4 appears in the am- 
plitude with a new interaction point, which may happen before or after the 
previous interaction point, thus the points of interaction lie chaotically on 
the time axis. However, from the point of view of a stationary process, which 
creates the system with a given quantized energy state in the limit of long 
time interval, one may think of an averaged progressive time and averaged 
trajectories of constituents, where stochastic time fluctuations are dealt with 
in a well defined averaging process. In this way the time-fluctuating relativis- 
tic trajectories are averaged into stationary time-ordered trajectories, similar 
to the quantum mechanical ones, where fluctuations are allowed for spacial 
coordinates. Correspondingly one can write 

z 4 (r) = z 4 (r) + z 4 (r), (5) 

where 24 (r) = £e = 2ur is the averaged time, proportional to the proper 
time, while the fluctuating time 2:4 (r) can be written as a sum of one-step 
fluctuations: 

n N 

z 4 ( T ) = Az 4 (fc), t = ne, Ne = s, ^ Az 4 {k) = 0. (6) 

k=l k=l 

The proper time s is expressed via the total Euclidean time T = x 4 — y± and 
the new variable u, 

s = T/2u>, (7) 
and hence the scalar Green's function ([T]) can be rewritten in the form 

duj 



roc W/,i 

g(x,y) = T / D 3 ze- K ^m Xl y)} Az4 . (8) 
J 2uj z 



Here 



T , I lo m 2 00 ( dz ^ 



*M = y o + _ + j, (9 ) 

while 

(Q(x, y))Az 4 = -DAz 4 exp[ze J A i (z(tE),tE+Z4 : )dzi+ie J A^dtE+ie J AidAz^]. 

(10) 
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DAz * = / it- II rrzz ex P 1 L 



-fc=i 



2tt 

The result of integration in ( TrTI) can be written as 
(*(z,2/))a, 4 



ip^Az^k) - ?±-LL- + zeAz 4 (fc)A 4 



2ttT 



(12) 



where y) is the averaged Wilson line, augmented by the fluctuation, 



$(x, y) = exp 



/•X rX4 

ie Ai(z(t E ),t E )dzi + ie AJz(t E ),t E )dt E 

Jy h*. 



exp(AS). 
(13) 

In the simplest case of the free scalar Green's function = and 
$(x, y) = 1, hence 



/ \ T f°° dbj 3 , 
0o(z,y) = \h~ — 7=(- D ^)xye 

V o7T JO WVW 



-*(«) 



8tt 2 T 



dcuexp 



m 2 T x-y 2 wT 

i U-u 

2u 2T 2 



n-^i(H, « 2 = T 2 +(x-y) 2 , 
(14) 



where i^i (x) is the Bessel function of the second kind. 

At this point the role of u becomes clear, since for large T the integral in 
(fl4|) can be taken by the stationary point method with the action 



s(uT) rn*T (x-y) 2 dS(u,T) 



o. 



(15) 



U}=LUQ 



and 



mT 



OJQ 



x 2 - y) 2 + T 2 



— > m, for T ^> |x — y| 



and one finally obtains the standard answer for large T yielding the asymp- 
totics of the r.h.s. of (1141) at large T. 



9o(x,y) 



m 



47r 3/2 T 3/2 



exp(— mT). 



(16) 
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Figure 1: The time-fluctuating trajectory in the (-21,-24) plane. The points 
( z i(tE), Ie) are marked by circles and are connected by average trajectory, 
depicted by dotted line. 



Another form exploits the Hamiltonian in (Q3J), namely, one can use the 
relation 

J (D 3 z)xye- K ^ = (x|e-^ T |y> (17) 

where 

H(u J ) = p2 + m2 (18) 

Applying the stationary point method to the integrals ({14]) . ( IT7|) . one 
obtains at large T the energy eigenvalue 

dH(co) 



duj 



u = uj = 0; uj = \Jp 2 + m 2 . (19) 



From (1191) one can understand, that oj plays the role of a virtual particle 
energy, and the condition ffl9|) has the meaning of the energy shell condition. 
This interpretation holds also for the case of N particles with interaction, 
when the integrals aver n£Li d^i are involved. Note, that in this way Ui are 
not any more approximate einbein variables, as in our previous works see e.g. 



We now turn to the case of 7^ and remark, that A^(z } za) are func- 
tions of coordinates, which will be used later in the process of vacuum aver- 
aging, yielding points of interaction, correlators etc., but at this moment in 
( TT3l) . $(a;, y) is a set of all possible Wilson lines, obtained by time fluctuations 
with the weight, given in ( TTTT) . see Fig.l as an illustration. 

In one particular case, when z(tg) is fixed, i.e. the trajectory is parallel to 
the Zi axis, all Az± fluctuations are washed out, since all fluctuations cancel 
each other, 

exp(ie J Ai(z,Z4)dz4) = exp(ie J A^z, tE^dtg). (20) 
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The same would happen in the case of QCD, where again overlapping 
pieces of Wilson line cancel each other. We shall come back to the problem of 
fluctuating Wilson lines, when we consider gauge invariant two-body Green's 
functions. 

In the case of the white system of the quark and antiquark of opposite 
charges, one must start with the one-body Green's function 



Sfay) = (m.+d-igA-ie^-i ee (rm+D^ = (mi-D®){Tr$-(D®)>)£. 

(21) 

The path-integral representation for Si is [5] 



S, 
where 



POO 

(x,y) = K - £>«) / d Si (Dz) xy e- Ki $®(x,y) = (m t - D®)Gi(x,y), 

J 

(22) 



1 rsi (dz^ X 2 



= + i /„ "Hit • < 23 > 



*<? (x,y)=P A P F exp (ig f A„dzf + ie, f A\ 

V Jy Jy 



X 



exp (J dna^igF^ + etB^)) . (24) 

Here F^ u and B^ v are correspondingly gluon and cm. field tensors, Pa, Pf 
are ordering operators, = ^(7^ — 7^7^)- Eqs- (12 lti24p hold for the 
quark, i = 1, while for the antiquark one should reverse the signs of and 
g. In explicit form one writes 

( <rH <tE \ / <xB \ 

<W> ~~ ^ c-E crH J ' ~ ^ aB J' [Zb) 

The two-body gig 2 Green's function can be written as [TTJ l2"Tj 

/•OO /-OO 

G qiij2 (x,y)= ds x ds 2 (Dz^) xy (Dz {2) ) xy (TW a (A))AX 



x exp(zei /* A^dzJ _ ^ 2 f Afdzf + ei dn(o-B) - e 2 dr 2 (<rB)), 

(26) 

where 

T = tr(ri(mi-Di)r 2 (m2-D 2 )), (27) 
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"tr" is the trace over Dirac and color indices acting on all terms. Here 
(W a (A)) is the closed Wilson loop with the spin insertions and one should 
have in mind, that color and e.m. spin insertions in general do not commute, 
which should be taken into account when computing spin-dependent part of 
interaction, see [36J, in f )26|) this fact was disregarded. 

fS 2 



W a (A) = P a P F exp 



, (28) 

It is important, that the physically meaningful result for the Green's 
function is obtained by two different averaging procedures applied to the 
total Wilson loop W = $^(x,y)$£\y,x) : 

1) one should overage W over all time fluctuations ; 

2) one should average W over all nonperturbative (np) and perturbative 
(pert) field configurations with the weight, given by the standard QED+QCD 
field actions ,so that the final result is 

((W)) = ((W) Az4 ) AAM . (29) 

However, the class of processes of interest in QCD is very wide, since any 
process, starting and finishing with definite hadron states, such as formfac- 
tors, decays, hadron reactions, needs an explicit definition of initial and final 
states as eigenstates of the Hamiltonian H(u\, 0J2), and therefore can use the 
formalism, discussed in this paper. 

It is clear, that in the fluctuation averaging (W)^ Zi the result is an average 
Wilson loop, passing through the points {z(t E )+Az(t E ), t E +At E }, t E e(0, T), 
where Az(t E ), At E depend on T, mi, m 2 and also on the concrete field config- 
uration, which will be averaged in the next averaging process, over vacuum 
fields. 

One can estimate the average time fluctuation At E in the case of the free 
relativistic particle propagation. 

E.g. assuming the correlation function to have the form 

- 4 2) ) = exp ( JZ \~f r ) , (30) 

and taking into account time fluctuations = t E ^ +z 4 (t E ) } and integrating 
over Az 4 , one obtains the increase of the correlation time 

(A^-) 2 (Az)L 4 = (Az) 2 + £ ~ (A,-) 2 + |-. (31) 
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However this result is an artefact of the not accurate definition of the 
path-integration measure, when at the ends of the time interval Atg the 
path can change the direction, implying infinite time derivative. Imposing 
a proper condition on the magnitude of the derivative, i.e. with smooth 
trajectories, the result would be different. From the point of view of the 
relation AM At > 1, one can in principle calculate however accurate values 
of masses M for large T, and only the coupling to decay channels, i.e. the 
width T should put a lower limit on the accuracy AM. 

It is interesting, how this problem occurs in our path-integral formalism. 
Indeed, the basic dynamics which is contained in {(W)), when time fluctua- 
tion is supported by interaction, Eq. ( 13"U|) . can be described by the diagram 
in Fig. 2. Now, from the point of view of Hamiltonian dynamics with the 
trace of the hypersurface, shown in Fig. 2 by a dotted line, the Hamiltonian 
becomes a matrix, with Fock states, numerating columns and rows, 



where nondiagonal elements are transition operators and diagonal ones define 
dynamics (and masses) of more and more complex systems. Therefore, e.g. 
V13 is responsible for the decay qq — > (qq)\ + (<?<?) 2, and hence defines the 
accuracy of possible mass determination of the incident state (qq). 

The wave functions of (132 p are actually Fock columns of different states, 
e.g. i^qq, ^(qq)(qq),...}, and therefore the (qq) eigenstates {$q^} are not any 
more an orthonormal set of states. As we shall see, the eigenstates ty n (ui, u 2 ) 
will not be orthonormal on the energy shells (u)i, 00^ ) , different for each n. 
In this way going from 4d path integral to the relativistic Hamiltonian 3d for- 
malism one naturally meets the many-channel Hamiltonian, where diagonal 
elements correspond to the fluctuation-averaged trajectories. 

3 From path integral to instantaneous dy- 
namics 

As a result of two averaging processes; time fluctuation and vacuum averaging 
the basic dynamical input of the resulting 3d path integral - the doubly 




(32) 
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I 



II 




Figure 2: The vacuum averaged Wilson lines, displaying pair creation in the 
time fluctuation process. The hypersurface traces I and II mark the double 
quark pair state of the Hamiltonian 



renormalized Wilson loop can be written as 

((W)) = Z w ex P ( [ -^J J d7T^(l)d7T Xa (2) [g 2 (F^(l)F Xa (2)) + 



e 2 (F$(l)Fg(2))]+0(FFF)} 



(33) 



where 



dt (1) dt {2) 

(1) al E _ (2) al E 



ft i' 



2uj? 



and the integration ds^is done over the minimal area S m i n inside the time- 
averaged trajectories of quark and antiquark L\ and L2. Note, that in addi- 
tion to the time-fluctuation smearing discussed above, there is also nonper- 
turbative smearing, provided by the np field correlators. 

Indeed, the quadratic (Gaussian) color field correlators can be written as 

rna 



n 2 I f)D E \ 

2 r «TrE t (x)*E J (y)&)) = \D E (u) + D? (u) + u^j 



+ UiUj 



dDf 
du 2 
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n 2 ( f)D H \ 8D H 

j r «TrH i (x)*H j (y)&)) = % [D h (u) + D»{u) + u 2 ^J - u iUj -^, 

2 8D EH 
^((TrH^Ejiy)^}) = e ijk u A u k ^^, (34) 

where D E , D H are purely np correlators, and D E ' H contain perturbative part. 
The same type of equations, but with replacement 4 — > e 2 and keeping only 

E H 

D 1 ' holds also for e.m. correlators. Note, that at zero temperature color- 
electric and colormagnetic correlators coincide, note also that np correlators 
D E , D H are due to Euclidean vacuum fields. 

The explicit form of perturbative correlators D E,H to lowest order in a s 



is 

D E (x) = D?(x) = ^ + 0(a 2 ), (35) 

while for e.m. fields one should replace |a s —> a. 

At this point it is important to realize, that the correlators depend on 
space and time intervals, e.g. D(z^ — z^ 2 \t^ — ^e) and ((W)} in Eq. 
(|34p even after fluctuation averaging implies nonlocal in time dynamics, e.g. 
the term / / e z {F^J{l)F^(2))ds^ u {l)ds\ a {2) stands actually for a photon 
exchange diagram. We are now going to replace this time nonlocal interaction 
by the instantaneous one, which is easily done in the correlator language, 
simply by integrating in fl34|) all correlators over time differences, — , 



dtfi dtfi = dt £jd(t £ — /. j,- ), dtjr; = d- 



f(l)j + (2) _ AtJuQ) _ + W\ - + f E 



It is important, that the main part of our interaction, the confining in- 
teraction, is ensured by the correlator D E (\/t 2 + r 2 ), which has a very small 
correlation length A, as was shown on the lattice [36] and analytically [37J. 
D E {t) ~ e~*/ A ,t > A, A ~ 0.1 fm, and therefore the transition to the in- 
stantaneous dynamics is done on small averaging interval At ~ A. Therefore 
for all processes with momentum (energy) transfer AQ, satisfying AQX < 1, 
this transition of np confining mechanism to the instantaneous dynamics is 
allowable. The case of gluon exchange is similar to the Coulomb interac- 
tion, where the instantaneous approximation in the Bethe-Salpeter equation 
is known as the Salpeter equation and is widely used in the literature. We 
shall mostly use the one-gluon exchange (OGE) interaction as a perturbation 
and therefore our transition to the instantaneous dynamics is justified. 
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For the case of the zero angular momentum (see in [23] the general deriva- 
tion) one can write for the instantaneous straight line w^(t,/3) = zjp(t)j3 = 
z^(t)(l -0),O<0<1, and e.g. ds^ = (z^(t) - z^(t))d/3dt. 

For zero angular momentum one can simplify the integration over the 
area of the minimal surface in fl33|) and obtain the result, neglecting spin- 
containing terms in (133]) for the moment, 

((W0> = Z w eM~ [ T [V (r(t E ))])dt E ), 
Jo 

where r(t E ) = \z-i(t E ) - z 2 (t E )\, and 

V (r) = V conf (r) + V OG E(r), 

V C onf(r) = 2r dX dvD{\, v) -> or, (r -»■ oo), 
Jo Jo 

/*oo /*oo 

a = 2 dv d\D(u,X), 
Jo Jo 

Voce = f XdX f°° duD^X, v) = -\*L 
Jo Jo or 

As a result one can write for the product of qq Green's functions (we 

omit renormalization Z factors, Fock amplitude coefficients, and ordering 

operators for simplicity) 



(36) 

(37) 
(38) 

(39) 
(40) 




(41) 

where A = Kx(ui) + K 2 (u 2 ) + J V (r(t E ))dt E , and 



We can also introduce here the two-body 3d Hamiltonian H(u>i, uj 2 , Pi, p 2 ) 
and rewrite ( l4"Tj) as 



{m\- Dl)(ml- Dl) 



j-y 



T 



57T JO 



duii 

3/2 
"1 



dou 2 

3/2 



x e 



-H(wi Wl ,p 1 ,p 2 )T| \ 



(42) 
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where H is obtained in a standard way from the action A(u x ,U2, z 1; z 2 ) (we 
omit all e.m. fields except for external magnetic fields B) 

2 (p(0 _ &(b x z«)) 2 + mf + uf - e^B 
# = £- -i +V Q (r)+V ss +AM SE (43) 

i=l ZUJ i 

and Vq is given in ( 137|) . The spin- dependent part of H, V ss is obtained per- 
turbatively from a^F^ terms in ( 128|) . and is calculated in the presence of 
m.f . in [36] . It is considered as a perturbative correction and is a relativistic 
generalization of the standard hyperfine interaction, 

V ss (r) = — ^ / (of v F, u {x)o$F pX {y))d{x, - y A ). 

7« FfoUM i 



Its explicit form is given in [38] . Finally, the correction — — F< £}^ EMl ; where 
i refers to the same quark (antiquark) yields the spin-independent self-energy 
correction AM$e which was calculated earlier [21] and for zero mass quarks 
and no m.f. is 

a „ r 3(7 3<T . . 

AM SE = — . (44) 

Z'KUJl 27TL02 

For the case of nonzero m.f. the resulting AM$e is given in [38]. We can now 
write the total Green's function of gi<j2 system, denoting by Y the product 
of projection operators Y = T(rrii — Di)T(m2 — D2), 

mi — -Di = mi -zpi = mi+a;i74 — zp7, m 2 - D 2 = m 2 ~ ^274 - «P7, (45) 

where p is the quark 3 momentum in the cm. system. 
As a result one has 



y>(x-y)G(x,?/)= y>(x-y)tr 



4F r 



{m\-Dl){ml-D 



2\ 



XI.J 



-|f%f%W(x|e-(— P,P^| y >, (46) 

Z7T JO w x ' Jo UJ 2 

We have used in fT46|) the relations 4(F) = tr(r(mi — if>i)T{rri2 — ^3 2 )), 
and neglect spin dependent terms in H; we have taken into account, that _D M 
acting on Wilson line, i.e. exp(ig J x A^dz^A yields exp(ig J x A^dz^d^A. 
The cm. projection of the Green's function yields 

[ rf 3 (x-y)(x|e-^ 1 ^'P 1 'P^ T | y ) =Y,<pl{0)e- Mn(fJl *» )T , (47) 
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see Appendix 1 for explicit separation of relative coordinates, Eq. (1A.1 9ft - 
(1A.1 12|) . Here M n {uj\,uj2) is the eigenvalue of H(uj\, cj 2 , Pi, P2) i n the cm. 
system, where P = p l + p 2 = 0; p 1 = p = — p 2 . 

The integrals over du\, du 2 for T — > 00 can be performed by the stationary 
point method, namely one has 



o; x Jo w 2 « 



^(0)(K) 



^4 0) 4 V(4 0) m;(i))(4 0) m;(2)) 



where 



0, Mji) 



<9M n (wi,u; 2 ) 



,(°) 



,(°) 



(48) 



(49) 



and we have neglected the mixed terms q^™ 2 for simplicity, however should 
keep them in concrete calculations: see exact result in Appendix 1. Com- 
paring the results f|4T|) . (jUJ) with the definitions of quark decay constants 

fn 
JTi 



G r (x)rf 3 x = E/ d 3 x(0\j r \n)(n\ Jr \0)e' 



Px-M n T 



d 3 P 



2M n (2Tif 



E £ r <S> e r 



(Mnf£ 

2M„ 



n\2 



-M n T 



where for T = 7^, 7^75 



fc=l,2,3 



(50) 



(51) 



and e r = 1 for T = 1,75, one obtains the expression for fp (to lowest order 
in V ss ) 



n\2 



^''r;' ," 2 . «. s v / (-S 0, m;;(i))(4 0, m;(2)), (52) 



, , (0) , , (0) M £ 



It is interesting, that numerical estimates using ( 152]) and (IA.1 20 p are 
close to those, obtained in 



14 



4 Relativistic Hamiltonians of a meson in mag- 
netic field 

The resulting relativistic Hamiltonian in the instantaneous limit ia given in 
and can be written as 



2 ( (Q _ a (B x zW)) 2 + rnj + uj - e^B (1) (2) (1) ( 

H-}_^ " 7, ^ U ( z z ; j ; ,Wi,W2) 

i=l ^ 

(53) 

where 

C/ = y ( r ) + V aa + AM SE (54) 

We shall be interested in the spectrum of the q\q 2 system in the magnetic 
field B, but before that we shall test the general form of the Hamiltonian 
H(oji,U2,PiP 2 ) and its eigenvalues, obtained at the stationary point values 

We start with the case of B = and U = — i^i^^pyr - Separating the 
total and relative momenta and coordinates, 

R= " 1Z(1)+ " 2Z(2) , r, = *V-*V' t W = IA, (55) 
LOi + Ul2 i or) 

and P = p^ 1 ) + p^ 2 -*, one obtains in () with B = 

H= ? ^ - + $ + U{ r)) + £ !±t«L. (56) 
2( Wl +w 2 ) 2u K,J ^ 2ui y ! 

1. The first example is the relativistic electron with mass m\ in the Coulomb 
field of a heavy atom of mass m 2 with charge Ze, ^(77) = — — . For 
P = one has for the ground state 

Minimizing in oj\ for m 2 3> mi one obtains 



M ^m 2 + m 1 ^l- (Za) 2 , (58) 
which coincides with the exact answer from the Dirac equation. 
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2. As a second example we consider electron-positron system, then from 
the same Hamiltonian ( 15 6 p for P = and m\ = m 2 = m one obtains 
after minimization 



M = 2my 1 - — « 2m — , (59) 

which looks correct, at least in the expansion in a. 

3. In the next example we consider the noninteracting qiq2 system in 
constant magnetic field B along the z axis. For U = one can solve 
one-body problem for each quark in m.f. with the result for the lowest 
Landau levels (LLL) 



M(u) 1 ,w 2 ) = 2^ ^— (60) 

i 1 

and after minimization one has 



M(4 0) ,4 0) ) = £ V(P^) 2 + ml + eB{2m + 1) - e^B (61) 



We turn now to the general case of the qiq2 system and consider first 
the case of a neutral system, e\ = — e 2 = e. In terms of total and relative 
momenta the Hamiltonian has the form 



(62) 



H 



B 



+ 



2wi 
1 



2coo 



^P + 7 r-^Bx(R+^r7) 

LO2 2 oji 



LO\ 2 LO2 



+ 



i=l,2 



2CJ,; 



(63) 
(64) 
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The R dependence for (I6"2j) in the case, when e\ = — e 2 can be factorized 
out in the way, discovered long ago 



tf(R, rj) = tp(rj) exp(zPR --(Bx ry)R) (65) 
and for (p(rj) one obtains the equation: 

(Hq + Hv + UMv) = M<p(ri), (66) 

where H is 

^^(-^ + T< Bx ") 2 )' (67) 

In this way Hq adds to the confining well with OGE and other terms also 
an oscillator potential. 

One can replace for simplicity the linear confining term by the oscillator 
potential, V con f = err] — >■ V con f = | (^+7)> where 7 satisfies stationary 
point condition ^| 7=70 = 0, which ensures some 5% accuracy of this re- 
placement. Then the lowest eigevalue M of the basic part of Hamiltonian, 
H — H + H a + V con f, is 

Mfa, o; 2 , 7 ) = e B±iB . + £ mf + ^~ e ^ )B ) (68) 

where e\ = e = — e 2 , and 



7<7 

+ -y, (69) 

We turn now to the case of charged two-body system in m.f., and here 
one can consider two different situation. In the first case, when e\ = e 2 = e 
and also mi = m 2 (and hence = o; 2 ) one can do an exact factorization 
of R and 77. 



"nj_,n z 



2Q 



/e 2 5 2 + —(2n ± + 1) + y — + 3 



P 2 eP(BxR) e 2 „ N2 tt 2 e 2 (B x r/) 2 (o; 3 + o; 3 ) A rr , 

Fb = of I S 4: -+^(BxR) 2 +— + 1 " 11 ^+Aff B ( Wl , 

2(o>i + o; 2j ) o>i + o; 2 80; 2o> o\u\ + o; 2j ro;io; 2 

(70) 

i^o- is given in and AHg is 
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AffaKo*) = . f,r(B x „) - ^^(B x R)- 

Wi^l^i + u; 2 ) 2 UJ1UJ2 2 

- ^^P(B x rt) + ? ^ 2 - ^i 2 ) g (B x R)(B x } (n) 
For cji =6^2) Ai?e vanishes and the Hamiltonian has the form 



TT P 2 e(P(B x R)) e 2 ,„ _ tt 2 e 2 /TI x2 
H — 11 + — B x R 2 + — + — — (B x r; 2 + 

Aui 2co Alo u 16a; 

2m 2 + 2w 2 - e(o-i + cr 2 )B <r /n 2 \ Tr 
+ — +2 I 7 +7 j + ^oge + V ss + AM 5S . (72) 

The lowest eigenvalues of the Hamiltonian ( 1721 ) are 

777 2 4- L} 2 

M = — - + (V ge) + (V ss ) + (AM SE } + 



eB , nnT feB\ 2 2a . , 1. [2a e{a x + cr 2 )B 7 (7 

+— 2iV ± +l + X hi- + 2n ± +l +(n„+- W 1 1 2; 

2a; \ \2o;/ 70W 2 y 7 o; 2a; 2 

(73) 

We now turn to the general case of a charged gig 2 system, when e± ^ 
e 2 , and write the full instantaneous Hamiltonian as in (1621) - (1641) . but with 
arbitrary e\ and e 2 , ei + e 2 = e and e is the total charge of the meson. 

In this case the simple factorization form fl65l) does not work, and one 
must instead to make a first step towards factorization, namely one must 
associate the cm. motion in m.f. with the total charge e of the system. This 
is done in the following form, discussed previously in 



*(77,R) = exp(TX77,R), (74) 



r = PR-|(Bxr ? )R, e= 6 -^ (75) 



and the resulting Hamiltonian from the relation H Q ty = exp(iT)H' (p, is 
2(w!+w 2 ) 2 2a) 2 F 
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+X Lri BL v + X 1 P(B xrj)+ X 2 (B x R) • (B x 77) + 



+ X 3 ir(B x R) + 



m\ + uj\ ml + u 2 



n 2 R = B 



n 



B* 



2wi 
2 ( ei + e 2 ) : 

2 



2w 9 



(eio;2 + ecui 



(e 2 u>i - ew 2 ) 



2wi 2w 2 

Here all coefficients X^i = 1,2, 3) given explicitly in the Appendix 2 of 
Treating the terms Xi,X 2 ,X 3 as a perturbation AMx, 



(76) 

(77) 
(78) 



AM X = (XxP{B x 77) + X 2 (B xR)(Bxjj) + X 3 7r(B x R)), (79) 

one can write the total energy eigenvalues of the Hamiltonian H in (1681) 

as 

M W = M (o) + M (o) ^ + AMx + ^ ( g0 ) 

where 2 

M {0) (P) = 7^-^ r + n R (2n Rx + 1) + X LP L P B, (81) 

(77) is the eigenvalue of the operator H n , 

tt 2 unw ± 

= — H h X Lri Bl, v + V co „/ + V ge- (82) 

We have written above the most general forms of instantaneous Hamilto- 
nians in the external m.f. It is seen, that to a good accuracy the dynamical 
contributions of e.m. and color fields can be separated, except in the OGE 
and spin-dependent terms, and as shown in [35] , the m.f. contribution to the 
both terms is decisive at large eB. 
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5 Discussion of results 



We have started with the general 4d proper-time path integral for the Green's 
function of a quark and an antiquark in gluonic (A^, F^ u ) and e.m. (A*, B^ u ) 
fields. These fields are contained in the generalized Wilson loop W with 
inclusion of spin-field operators (a IJlll (F f j iV + B^)). 

After vacuum averaging procedure in the partition function, the averaged 
Wilson loop (W) A ^ A (e) contains all possible interactions, including internal 
quark loops from the terms tr ln(m 2 — Df) in the partition function. 

As a first step we have traded the particle proper times for the Euclidean 
(ordering) times t$ , and performed path integration over fourth particle 
coordinates £4,^4, which is physically the time fluctuations around t^\t^. 
We have shown, that this time-fluctuation integration leads to the 3d path 
integrals with the action (or Hamiltonian in the Hamiltonian form of path 
integral) which is a matrix in the Fock states. The resulting 3d path integrals 
are integrals over new parameters u±,u 2 , and the spectrum of the qiq 2 system 
can be found for large times by a stationary point procedure in Ui,u 2 . 

In this way one is going from the 4d formalism to the multichannel 3d 
formalism with an additional u- integration for each particle. 

As a next step we have observed that the interaction appearing in the av- 
eraged Wilson loop, (W) AA (e), has the form of field correlators (F fll/ (x)F\ a (y)) , 
(B^^Bxaiy)), and the first correlator has a very small correlation length 
A ~ 0.1 fm (found on the lattice [SB] and in analytic calculations [37]). This 
allows to go over to the instantaneous dynamics, when the bolocal or mul- 
tilocal) interaction (F(x)F(y)) is replaced by the time-averaged potentials 
V^(x — y) = / d(x4 — y±)(F(x)F(y)) , and this is valid when the basic pa- 
rameter, defining the quark trajectory, string tension cr satisfies crA 2 <C 1, so 
that typical time length on trajectory t ~ 775 is much larger than A. Note, 
that this condition is opposite to the one, used for validity of the OPE and 
QCD sum rules. 

As a result one obtains the instantaneous relativistic Hamiltonian H(ui, u 2 ) 
depending on two parameter u\,u 2 (for the (<Zi?2) Hamiltonian matrix ele- 
ment) and the actual spectrum is obtained from the eigenvalues M n (u; 1 , u 2 ) 
at the stationary points lo^ , uj^ . Note, that these points are different for 
different n = 0, 1, 2... 

We have checked the results in section 4 for several simple systems and 
found good agreement with known results. Moreover, this formalism for 
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eigenvalues was being used for more than 20 years in many papers, a small 
part of which was cited here, and the results in all systems, mesons, baryons, 
hybrids and glueballs are well compared with experimental and lattice ones. 

The important new element in this paper is the rigorous derivation of the 
integral representation for the (gi^) Green's function Eqs. ( HT|) . f )42|) . (jlSjl . 
which gives a new meaning to the parameters U\,U2, and allows not only 
calculate spectrum, but also the Green's function itself. 

As an important application of the developed formalism, we have derived 
in section 4 the explicit form of Hamiltonians of the (#1(72) system in the 
constant m.f. B, and defined the main part of the spectrum for neutral and 
charged mesons. 

These results have been used for the explicit numerical evaluation of the 
p-meson spectra in [JT] , which are in reasonable agreement with existing 
lattice data. Moreover, the same formalism was extensively exploited in [35] 
for calculation of chiral condensate, and in [33] for magnetic moments. 

Actually, the field of possible applications of our method in QCD and 
QED is enormous, and the method is especially simple in the cases, when only 
spectral properties are of interest. This is clearly seen, when one compares 
this method with the Bethe-Salpeter equation. In the last case one is facing 
the problems of the relative time and insufficiency of the ladder kernel already 
in the QED case. 

In the QCD case the use of the Bethe-Salpeter equation is in addition 
associated with the vector propagator form of confinement, which is physi- 
cally not consistent, or with some phenomeno logical form, and in this way 
the method loses its fundamental character. On the contrary, the very short- 
correlation property of confinement suits perfectly to establish the validity of 
instantaneous Hamiltonian formalism and allows for an accurate and simple 
procedure. 

The author is grateful to M.A.Andreichikov, A.M.Badalian, S.I.Godunov, 
B.O.Kerbikov, V.D.Orlovsky, A.E.Shabad, M.I.Vysotsky for many useful dis- 
cussions. 



Appendix 1. 

Derivation of the general expression for the qiqz Green's function 
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We start with the general definition for the q\q 2 Green's function in the 
vacuum gluonic and external e.m. fields 

G(x,y) = (trTS 1 (x,y)TS 2 (y,x)) A = 

= <^^F#r<^4^ = (A.i i) 

m i " D i (m| - D 2 ) 

roo roc 

= 4/ d Sl / ^(Z) 4 ^ 1 ^ 4 ^ 2 ^^^ 1 ^^^}, (A.l 2) 

JO JO ' 

where (I^Wf) = jtr[T(m 1 - zpx)f (m 2 - ip 2 ) (W F ) A ], and W> = ({W}}, given 
in (1551) ; the spin operator ordering in (IA.1 2\i is not written explicitly. Ne- 
glecting spin dependence, one has a purely scalar function W F , which is 
proportional to a unit (4 x 4) matrix. 

Introducing the effective energies u>i = T = |x 4 — y 4 |, one can rewrite 
(lAXTD as 

G(x,y) = L r%^(^^ J D 3 ^)xye-^ ( - ) -^ ( - ) (W))A Z4 , 
2vr Jo uj-l' 2 co^ 

(A.l 3) 

and we have taken into account, that 



(Dz^Dzi% 4m (YW F )e " u V - / ^ v ' * J = ^^((YW F )) Azi 

(A.l 4) 

Here ((YW F )) & Z4 corresponds to the time-fluctuating Wilson loop aver- 
age, as in Fig. 1, renormalized and normalized by the condition 

((YW F )) Az4 (g = e = 0) = l (A.l 5) 

We omit in what follows the Fock column structure of the corresponding 
particle contents in our averaged Wilson loop ((YW F ))fi with the corre- 
sponding Zi factors for each Fock line and concentrate on the simplest case 
of one renormalized closed {q\q 2 ) loop depending on t F , as shown in Fig. 2. 
In the neutral case e\ = — e 2 , ((YW F )) a 24 depends only on coordinate differ- 
ences r](t E ) = z^(t E ) — zP\t E ), defined at the same moment t E , and one 
can proceed integrating out the cm. motion. Ki,K 2 in ( 1A.1 3ft are 



K 1 (u 1 )+K 2 (u 2 )=[^- - + -^ 2 -)T+ dt E 

V lLO\ 2u 2 I Jo 



u x fdz^V u 2 u<zpy 

2 \ dt E J + 2 \ dt E 
(A.l 6) 
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Introducing now the coordinates 

rj{t E )=^-^\ p(t E ) = ^-zV>(t E ) + ^—*V>(t E ) (A.17) 



one can rewrite the last term in (1A.1 61) as 



and the path integral (D 3 D 3 z' 2 ))xy as 

fD 3 z( 1 )D s z( 2 )u= /" rf3pi /" rf3p2 P i Pi(E Az(1) -( x -y))+^p 2 (E Az(2) -( x -y))> 

1 ;xy y (2tt) 3 y (2vr) 3 

d 3 Az^ d 3 Az^ 

X (i^P (4^ = (^xyPUo, (A.1 9) 

Where 

c^WS?^^^, (A , 10) 

In absence of external magnetic field, which acts on cm. coordinate p, 
it is convenient to consider the P = projection of the Green's function 

/ G(x,y)d 3 (x-y) = L f % r^(D 3 rj) 00 e- K ^((YW F )) Az , = 

J Z7T JO iiU JO (.)_' 



where 



A"" = + + A ** 2 bs ■ (A- 1 13) 
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(0| e -" T |0> = J2 |^n(0)| 2 e- M " (wi ' W2)T . (A.l 14) 



n=0 



Here <p n (0) = ip n (u)i, co 2 , ?7)|™ = o> andM n (wi,w 2 ) is the eigenvalue of the Hamil- 
tonian 

H = H(ui,u 2 ), H^ n = M n (uj u uj 2 )(fn- (A.l 15) 
Assuming, that ((Wp))& Zi can be represented as 

((W F )) Az4 = exp(- J V(r),oj)dt E ), (A.l 16) 
the Hamiltonian can be written in the form 

u; 2 ) = £ + ^ + Vfo, w lf wa) (A.l 17) 

i=l 2 ^ W 



At this point one can define the so-called quark decay constants /r [30] . 

/ G(x,y)d 3 x-y) = Y,er<g>er Mni f )2 e~^ T (A.l 18) 

where £p = 1 for 5 and P channels, and £r = £^ for V, A channels, 

E ^? ) (?)4 fc) (ff) = ^-^ (A.i 19) 

fe=i,2,3 " 

and hence can be found from flA.l 1 21) as 



Here T on both sides is assumed to tend to oo, and one can calculate the 
integral on the r.h.s. of flA.l 20p by the stationary point method, 

,(o)> 



M n {u u u 2 ) = MM 0) M 0) ) + M^\u?\^) {U1 ~? )2 + 

(0) -+Mr ) (4 o) ,4 o) )(- 1 -^)(c 2 - 

(A.l 21) 



2 

+ iie^%f) ( ^° ))2 + m^A 0) )(^ - c{ 0) )(c 2 - 4 0) 
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where 



and 



M {ik) 



d 2 M n 
dujiduk 



(0) (0) 

L l J ft. 



0,i = 1,2. 



.(0) 



(A.l 22) 



(A.l 23) 



Doing the integration in (1A.1 20p with the help of ( 1A.1 21 [) one obtains 



(n) . 2 _ iV e 0O^(0) 



where 
with 



(c4 0) c4 0) )M n £ r , 



^ = \M 0) 4 0) ^n 



f2 n 



(«-/3) 2 + 7 2 
where we have denoted 

1 



a/3(a-/3) 2 ^[(a + Pf - 2(a - (3) 2 - 7 2 ] 



4[(a - /3) 2 + 7 2 ] 



(A.l 24) 



a 



7 



M (12) 



and finally 



M n = M n (uj[°\ lo { 2 0) ), Y = ^tr D (r(mi - ip 1 )T(m 2 - ip 2 ) 



(A.l 25) 



(A.l 26) 



and £t\d denotes trace over Dirac 4x4 indices. It is instructive to com- 
pare (I A.l 24[) with the old result, obtained in [30], using approximate path 
integrals over (DAu), 



{n)2 2iv c (i> 2 (o) 



(A.l 27) 



As one can see, comparing (1A.1 241) and (1A.1 27jl . in the first case (the 
time-fluctuation approach of the present paper) the factor , 1 should 



be equal to 2, for both expressions to coincide. In practice for the (gi^) state 
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made of zero mass quarks, mi = m 2 = 0, and with the total mass made of 
confining interaction, see [2S] for details, one has 

2 2 2 

M n {u u co 2 ) = £ ^i±^L + (2u;)- 1 /V 2 / 3 a n , a = 2, 338 (A.l 28) 

and for mi = m 2 = one obtains 

(fi wW 0) )- 1/2 = 3, (A.l 29) 

while for mi = 0, m 2 <C y^cr, the result is 

(fi c4 0) c4 0) )- 1/2 2.34. (A.l 30) 

This implies, that the quark decay constants /p obtained in the new 
method will be larger by (10-20)% as compared with previous calculations 
in jig. 
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